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The spatial distribution of fields and currents in confining theories can give direct 
evidence of dual superconductivity. We review the behavior of vortices in the lattice 



^ : Higgs effective theory We discuss the techniques for finding these properties and 

. calculating the superconductivity parameters in lattice simulations. We have seen 



dual Abrikosov vortices directly in pure U(l) and SU(2) and others have also seen 
them in SU(3). We review the duality transformation for U(l) in order to connect 
the U(l) results to a dual Higgs theory. In the non-Abelian cases the system appears 
D . to be near the borderline between type I and II. We also discuss the response of the 



persistent currents to external fields. 
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I. INTRODUCTION 



The physics of superconductivity has provided the inspiration for many of the ideas be- 
hind the color confinement mechanism in QCD. The most notable property is the behavior 
of a magnetic field imposed normal to the surface of such a material. For a type II supercon- 
ductor, the field drills holes in the superconductor forming a lattice of filaments of normal 
material defining the cores of Abrikosov vortices [Q. The longitudinal magnetic field is 
a maximum at the center and falls off exponentially in the direction transverse to the core 
over a distance denoted as the London penetration depth, A. One quantum of magnetic 
flux is carried by each vortex. Recent photographs using a scanning-tunneling-microscope 
10] and electron-holograph-interferometry @] show this flux confinement very dramatically. 

Dual superconductivity |Q is the analogous phenomenon in which persistent currents 
of magnetic monopoles form dual Abrikosov vortices which confine electric flux. Lattice 
gauge theory offers the prospect of exploring dual superconductivity in depth as a conflning 
mechanism. In four dimensional U(l) lattice gauge theory there is considerable indirect 
evidence from "bulk properties" of the vacuum such as the monopole density 0, monopole 
susceptibility 0, and static quark potential that Dirac magnetic monopoles are associated 
with the phenomenon of conflnement. 

The existence of a dual Abrikosov vortex between a static quark-antiquark pair leads to 
more direct evidence for conflnement. By studying the relationship between the monopole 
currents and flelds in the neighborhood of static sources we identifled the flux tube with a 
dual vortex in U(l) PJlO(|, and in SU(2) |ni|. Suzuki and collaborators reproduced the SU(2) 
results with better statistics and verifled the effect in SU(3) [jl2|. We have further studied 
the effect at flnite temperature and have seen the vortex disappear above the deconflning 
temperature [|13|. Cea and Cosmai |]14[ have measured the electric fleld proflles in a flux 
tube in U(l) and SU(2) identifying the shape with a dual vortex. 

In ordinary superconductivity, the primary issue is the spontaneous breaking of the 
electromagnetic U(l) gauge symmetry (SSB) signaled by the non-vanishing of the vacuum 
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expectation value of a charged field. An immediate consequence is that the curl of the vector 



potential is proportional to the curl of the electric current known as the London relation [15 



The London relation is violated only near the boundaries of superconducting material within 



a distance defined as the Ginzburg-Landau [16| coherence length, ^. Combining the London 



relation with Maxwell's equations gives mass to the electromagnetic field. The Meissner |T7 
effect and infinite conductivity follow. 

These relationships suggest numerous ways to search for signals of dual superconductiv- 
ity. Numerical methods may be successful in establishing some of these connections and not 
others. The focus of this work is to pursue the dual London equation which is central to the 
phenomenon. There are other promising approaches which we will also describe. 

The theory of superconductivity entails (i) the identification of what symmetry is broken 
and what are the relevant dynamical coordinates, (ii) the mechanism that leads to the 
the instability and hence SSB, e.g. BCS [|18| theory, and then (iii) an effective theory of 
the currents and fields in the broken phase, e.g. Ginzburg-Landau [|T6| [GL] theory, or 
equivalently the four dimensional generalization, the Higgs model [|I^] treated as an effective 
theory, which allows the calculation of the spatial consequences of the broken symmetry . 

Our goals are similar to these. We hope to identify relevant dynamical coordinates and 
a corresponding effective theory. Experience has shown that an analytic approach is very 
difficult. The duality transformation has not been achieved for the Wilson form of the action 
which we are using. But this transformation has been implemented in closed form in the 
closely related Villain |2y] form of the action by Polyakov [^] and Banks, Myerson and 



Kogut p2[ in the early development of this subject. 

Frolich and Marchetti [^] have developed these techniques further. They have identified 
the monopole field operator for the Villain action. The non-vanishing vacuum expectation 
value is identified as the order parameter for dual superconductivity. This has been studied 
numerically by PoUey and Wiese, and Polikarpov p^p5| . Ref are other closely related 

papers. 

Similar analytic methods have not been successful for the Wilson form for the U(l) action. 



However Del Debbio, Di Giacomo, Paffuti and Fieri p8| have identified the monopole field 
operator applicable for a more general U(l) action and found the vacuum expectation value 
to be a very strong signal for dual superconductivity. 

While the origin of confinement in U(l) lattice gauge theory is fairly clear, understanding 
confinement in non-Abelian SU(N) theories has been more difficult |2^. One promising 
approach is to fix the non-Abelian degrees of freedom in the maximal Abelian gauge as 
advocated by Schierholz and Suzuki and others [^,^], leaving a residual U(l)^~^ gauge 
freedom, with (A^-1) species of U(l) Dirac monopoles. Suzuki and collaborators p2| , ^ have 
laid much of the foundation for this lattice approach. The monopoles have been observed to 
be abundant in the confined phase and dilute in the (finite temperature) unconfined phase 
3g,|31],||,|3|g,|3|®|33®|39|. A number of contributions by Yee are described 



in a conference talk 40 



Del Debbio, Di Giacomo, Paffuti and Fieri p8| have applied their U(l) construction to 
the Abelian projected SU(2) and again saw a clear signal of dual superconductivity. They 
have considered a number of forms of the Abelian projection but have not yet implemented 
this method for the maximal Abelian gauge. 

Dual vortices have been seen for SU(2) as reported in Refs. [|TT],|T^|I3| . All these papers 
report the similar conclusion that the dual GL coherence length and the dual London 
penetration depth are roughly equal. For a type II dual superconductor = Xd/^d > 
1/v^, and a type I otherwise Maeda, Matsubara and Suzuki |^2| came to a similar 
conclusion for SU(3) from an entirely independent argument. For superconducting materials 
the difference between type I and type II is dramatic when a magnetic field is imposed. The 
penetration of magnetic fiux in a type I material leads to very complex patterns of normal and 
superconducting domains called an 'intermediate state structure'. Numerous photographs 



of this phenonenon are shown in a monograph by R. F. Huebener P3 |. 

Since U(l) lattice gauge theory does not have an non-trivial continuum limit, one must 
use the results with caution in describing confinement in QCD. However it is more than 



just an independent model of confinement. Yee has found using the demon method 



that for small (3 the Abelian projected action for SU(2) is essentially the U(l) action with 
Pu(i) = /^5C/{2)/2. Further this must and indeed Yee finds it does break down in an interesting 



way as j3 approaches the scaling region. See also Suzuki |32| for a different approach. 



In a long series of papers and a review article, Baker, Ball, Zachariasen |^ and co- 
workers have postulated a dual form of QCD in the continuum. Their formulation describes 
a non-Abelian dual superconductor and hence it confines color. They have calculated flux 
tubes, static quark potentials, temperature dependent effects and many other quantities in 
the tree approximation. Although their formulation does not focus on monopoles, there is 
a large potential overlap between their work and the lattice description developed here. 

In Section|I| we review the Abrikosov vortex in the Higgs theory |jl9[ which is the basis 
of the many interrelated signals of superconductivity. We will emphasize the lattice form of 
the Higgs model. The singularities giving vortices arise in an interesting way naturally on 
the lattice. 

In Sectionpn] we report the results of numerical studies of U(l) for correlations between 
a Wilson loop and a plaquette giving a very clear signal for the London penetration depth. 
These numerical results suggest that a dual form of the model should exhibit this phe- 
nomenon in a simple way. In Section|rV| we show that the simpler Villain ||2^ form of the 



action does allow an analytic treatment of the question and we review the results of Frolich 



and Marchetti p3[. Policy and Wiese [Q, and Polikarpov ||25|| . 

Results for SU(2) and SU(3) are given in Section^. The analysis involves an approximate 
treatment of the GL theory in order account for the behavior of the fields in the dual vortex. 

Fields and currents exist only near the boundaries of dual superconductors over a distance 
scale set by the London penetration depth A^. It is these spatially transient effects that 
provide spatial structures where the local properties of dual superconductivity can be studied 
and that is the focus of this work. One way to study these phenomena is to impose an 
external field. In Section|VI| we introduce an external field and discuss some of the prospects 
and problems. 

The lattice approach provides a very concrete way to visualize physical quantities as living 



on links, plaquettes, dual plaquettes etc. The use of the abstract language of differential 
forms is gaining use in the lattice literature. This formulation focus on k dimensional volumes 
and boundaries of the volumes in a systematic way. The technology in some ways removes 
one from the nitty gritty of the lattice but with some familiarity it actually enhances ones 



of this formulation. In the appendix we go one step further and give a concrete translation 
of differential forms into lattice difference operators. 



The Ginzburg-Landau effective theory of superconductivity is a time independent de- 
scription based on a complex order parameter, ip{x), where |'?/'(a;)p is the number density of 
superconducting charge carriers. In normal material = 0. We will focus on the special 

circumstance in which IV'I^?)! is constant in the superconducting phase and discontinuous at 
the surface. This is the extreme type II limit, also called the London limit. 

The Higgs |19| model is the four dimensional generalization of GL theory with the Higgs 
field corresponding to the GL order parameter ip{x). The London limit is obtained by 
constraining the Higgs field: 



In other words the London theory is equivalent to the effective non-linear Higgs theory. 



ability to see general relationships. Policy Weise give a brief introduction to the basics 



II. ABRIKOSOV VORTICES 




(1) 



A. London theory 



A concise statement of the London theory is contained in the relation 



A + \^J = 0; (V-1=0). 



(2) 
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If the charge density is zero, then in this gauge the electric field is given by —A and therefore 
E = X^J. This describes a perfect conductor and is just Newton's law for free carriers, 
eE = mv. By taking the curl of Eqn.(^ we obtain the condition for a perfect diamagnet. 

VxJ=-1r (3) 

This relation together with Ampere's law V x 5 = J gives V^B = B/\^ which implies that 
the magnetic field falls off in the interior of the superconductor with a skin depth A. 
Finally the fluxoid is given by the integral 

I {B + A^V X J)-nda= f (A + X^J) ■ dl = N— = Ne^ = (4) 
Js Jc e 

If the curve C is in a simply connected region of a superconductor, then = 0. However if 

the curve encircles a hole in the material then N need not be zero but must be an integer 

due to the quantization of flux. In an extreme type II Abrikosov vortex, a very small core is 

comprised of normal material. A single unit of magnetic flux of radius ~ A passes through 

the vortex. The fluxoid density B + A^V x J is zero everywhere except in the region of the 

normal material. In the limit in which the core is a delta function we obtain: 

5 + A^V X J = ^m5^{x±)n,. (5) 

Further if we use Ampere's law we can get an analytic expression for the B^ proflle of a 
vortex. 

B.{r^) - AVii?,K) = <^>,j\x^)- B, = l^KoirJX). (6) 

B. Lattice Higgs in the London limit 

I would like to review the salient features of the lattice Higgs theory in the tree approx- 
imation. It is interesting to see how the the London relation arises as a direct consequence 
of the spontaneous breaking of the U(l) gauge theory. Also we compare the simulations of 
pure gauge theories with this model interpreted as a dual theory. 
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Consider the lattice action: 

We need a background field in order to form an Abrikosov vortex. To keep as close as 
possible to our simulations we generate the field from a monopole loop, i.e. the dual analog 
of the Wilson loop. Hence we will digress to define the monopoles on the lattice using the 
DeGrand-Toussaint p construction. 

Consider the unit 3-volume on the lattice at fixed X4 as shown in Fig.|l]. The link angles 
are compact, —it < (f^ < it. The plaquette angle is also compact, —An < 6^1^ < Atc and 
defined 

ea?F^^{x) = d^y{x) = A+v9i.(x) - A+(p^(x), (8) 

where a is the lattice spacing and the difference operators are defined in the appendix. This 
measures the electromagnetic flux through the face. Consider a configuration in which all 
link angles on this cube are small compared to vr. Gauss' theorem applied to this cube then 
clearly gives zero total flux. Because of the 27r periodicity of the action we decompose the 
plaquette angle into two parts 

0^,u (x) = e^,u (x) + 27rn^^ (x) . (9) 

where — vr < 9^^, < n. If the four angles shown in Fig.|I] are adjusted so that e.g. 9^^, > vr 
then there is a discontinuous change in by — 27r and a compensating change in n^^. We 
can clearly choose the configuration that leaves the plaquette angles on all the other faces 
safely away from a discontinuity. We then define a Dirac string n^i, passing through this 
face (or better a Dirac sheet since the lattice is 4D). 9^i, measures the electromagnetic flux 
through the face. The important points are the following: 

• The volume shown contains a magnetic monopole as indicated by a net flux of 27r/e = 

Cm- 
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• Note that if one does a gauge transformation e.g. at the point x + e^'^) + e^'^^ then if 
one of the affected hnks picks up a discontinuity of 27r, the two contiguous plaquettes 
will be affected but the flux out of the volume is unchanged. The monopole is a gauge 
invariant construction. 

• This example created a monopole-antimonopole pair in the two space cubes separated 
by a lattice spacing in the direction. The same holds in the X4 direction (not shown 
in fig.|). 

• In the dual description, cubes in the original lattice are represented by links on the 
dual lattice orthogonal to the 3- volume. 

• Hence for n^i, = 1, (a Dirac sheet threading the plaquette) a closed 1x1 monopole 
loop has been created on the dual lattice in the x^^x^ plane forming the boundary of 
a Dirac sheet. 

• The monopole loops are closed, hence the current is conserved. 

This construction gives the following definition of the magnetic monopole current. 

— j;r(a^) = (10) 

Cm 

Which lives on the dual lattice. It satisfies the conservation law J™(x) = 0. 

We obtain a classical solution by minimizing the action. Our initial configuration con- 
tains a single closed magnetic monopole loop in the X3,X4 plane analogous to the Wilson 
loop projector in pure gauge simulations. This can be easily accomplished by choosing a 
configuration with a Dirac sheet threading all 1-2 plaquettes in a the 3-4 plane that spans 
the desired monopole loop. Our algorithm then rejects all updated links that change the 
initial monopole configurations. 

We accept configurations with Dirac sheets with no accompanying monopoles. A lump 
of Dirac sheets is rather easily created in a single local update. The update of a single 
link can produce six Dirac sheets in the plaquettes contiguous to this link. Each of these 

9 



produces a monopole loop on the dual lattice. One can see as shown in Fig.^ that the six 
dual monopole loops form the faces of a cube on the dual lattice in such a way that all the 
monopole currents on the edges of the cube cancel. 

We use the method of simulated annealing, slowly increasing jS, holding \^/a? = f3/K = 
constant, where A/a is the London penetration depth in lattice units. 

The electric current is given by 

— J'tix) = /m(0*(x)e*^'^(^V(a; + e^"^)). (11) 
en ^ 

Writing the Higgs field = p(x)e*'^*-^\ the curl of the current is hence given by 
(A+ J-(x) - A+ J-(a;)) _ 



BK a 

p{x)p{x + e^^^) sin[-cc'(a;) + ipf,{x) + uj{x + e^^^)] + 

p{x + e(^))p(x + e^'^) + e(")) sin[-uo{x + e^^^) + ^,{x + t^^^) + uo{x + e^^^ + e^^^^)] - 
p{x + e("))p(x + e(") + e^'^)) sin[-cj(x + e^")) + ^^{x + e^")) + uj{x + e^'^) + e^'^))] - 
p{x)p{x + e^'^)) sin[-c<;(x) + ^^(x) + uj{x + e^^^^)], (12) 

Compare this with the electromagnetic field tensor 

ea^F^^ = sm[ip^{x) + ip^{x + e^^^) - ip^{x + e^"^) - ip^{x)]. (13) 

If the U(l) gauge symmetry is spontaneously broken, these two quantities are equal which 
is the London relation. To be more precise: if (i) p is nonvanishing and independent of 
position (absorb the normalization into k) and (ii) 



thenB 



sin[^ + 2N7i] ^ e, (14) 



^ p a' i^pti^) - A+ J^(a:)) 27riV 1 1 

^ijiv = F^^ ^ = = A^e™ — , (15) 



^An apparent sign difference with Eqn.(y) is a notational one, see the appendix A2 
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T^vO? summed over the transverse plane is called the fluxoid and hence J^^i, is the fluxoid 
density. The magnetic charge is the Dirac monopole charge. Condition (i) is satisfied by 
freezing out the modulus of the Higgs field with the constraint Eqn.(|l]). and (ii) is satisfied 
if the variation of the field modulo 2'kN is much smaller than a lattice spacing. (Note that 
this singularity is not an isolated Dirac sheet. That would act differently, adding 27rn to 
both terms in Eqn. ([l5|) and giving a zero contribution on the right.) 

Eqn.(p!5l), together with Maxwell's equations, gives the lattice version of the Meissner 
effect: B — X^AB = 0, where = m~^; infinite conductivity E = X'^A^J (assuming p = 0); 
and an Abrikosov vortex 

B,-X\A^xJ\ = n'-^51^,. (16) 

Fig.0 shows the profile of the R.H.S. of eqn.([l^ and eqn. ([T3|) in the directions perpen- 
dicular to the 5x5 magnetic monopole loop on a 12^ lattice, with P/k = 1. We used 
the constrained form of the Higgs potential, eqn.(|l]), which corresponds to an extreme type 
II superconductor. The graphs show the expected behavior: i.e. the equality of the two 
quantities everywhere except at r = where they should differ by 2it. There are significant 
violations only at r = and smaller violations at r = 1 due to the breakdown of Eqn. (|I1|). 
We anticipate that the violations will decrease with increasing loop size. This solution is a 
possible way to treat the complications that are ignored in the continuum London theory 



described in Sec .III A 



Note that there is a sign change in the curlJ^ profile. Eqn.(^ shows that there must be a 
sign change since it must integrate to zero over the plane. The vortices have an exponential 
profile. All the current circulates about the core in the same sense. However since the 
profile falls faster than 1/r, the line integral around a small patch at the origin will have the 
opposite sign than the line integral around a small patch elsewhere. 

It is interesting to note that the interior surface spanned by the monopole loop is in 
the normal phase since the London relation is violated there. All other regions are super- 
conducting. This translates in the following sections in which the Wilson loop projects out 
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the normal phase in the plane spanned by the loop and the remaining regions are a dual 
superconductor. 

We chose the Higgs field \(t){x) \ = 1 for this discussion. If instead the field was subjected 
to a SSB potential 

Vh.,M{x)\') = Xmx)\' - v')\ (17) 

there are important differences in the solution. In the region near a superconducting normal 
boundary, one expects a soliton like deformation of the Higgs field, constrained to vanish in 
the normal region and favoring a non- vanishing expectation value deep inside the supercon- 
ducting side. This transition region defines the Ginzburg-Landau coherence length ^, which 
is the distance over which the order parameter rises to its asymptotic value. For distances 
into the superconducting material deeper than the coherence length, one expects the London 
relation to be satisfied. 



III. DUAL SUPERCONDUCTIVITY IN PURE U(l) GAUGE THEORY 

We now turn to the dynamical simulations in pure U(l) gauge theory given by 

S = (3 J2 (1-cos V(a;)). (18) 

X,fl>U 

We have given evidence ref. Q for the dual fiuxoid density relation signaling a dual Abrikosov 
vortex: 

a (X 

where *F^u = —^^fMuarFar- The magnetic monopole current is defined in Eqn.(|10|). 

The correlation of A x J™ with a Wilson loop gives this signal for the solenoidal behavior 
of the currents surrounding the electric fiux between oppositely charged particles. The 
correlators for the electric field E and x J"^ are 

< e^^^'e^^^ > < sin 9w sin 9p> 



< e^^w > < cos 6w > 



— (A~ X = < e''" ^i^- X Jr^t)^ > ^ <smew (A- x J,-), > 
Cm < e^^w- > < cos9w > 
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The notation (■ ■ ■)iat denotes an integer valued lattice quantity. 

The identification of sindp with the electric field is arbitrary, any quantity with the 
correct naive continuum limit will do. However Zach, Faber, Kainz and Skala made the 
interesting observation that this choice in conjuction with the Wilson action satisfies Gauss' 
law HBI . 



For out of the plane of the Wilson loop, the RHS of Eqn.(|l^) vanishes and the London 
relation is satisfied. Hence we interpret the Wilson loop as a projector onto the normal phase 
in the plane bounded by the loop and a dual superconductor elsewhere. 

Fig. ^ shows the lattice operators for the electric field and the curl of the magnetic 
monopole current. The longitudinal electric field, (a), is given by a t plaquette which is 
depicted by a bold line for fixed time. The curl of the magnetic monopole current, (b), is 
built from four 3-volumes which appear as squares since the time dimension is not shown. 
Passing through the center of each square is the link dual to the 3-volume. One takes the 
monopole number n in each 3-volume and associates the value ncm with the corresponding 
dual link. The 'line integral' around the dual plaquette completes the picture. Notice from 
this construction that E and A~ x J"^ take values at the same location within the unit cell 
of the lattice, both are indicated by the bold face line in the z direction. Both operators 
live on the same z, t plaquette. However x J™ is a hefty operator involving many links 
and as a consequence gives lattice artifacts for correlations at modest separations. 

Fig.^ shows the profile for the quantities making up the London relation on the same 
plot as for the effective lattice Higgs case, Fig.|^. We did 800 measurements on a 12^ lattice, 
P = 0.95, and a 3 x 3 Wilson loop was used to project onto the g, q sector using methods 
described in Ref. 0. 

Fig.^ shows the behavior of the dual fiuxoid density. We did a fit of the London 
relation, eqn.(^), using the R points and determined Xd = 0.49(3). The fit strongly 
favors the small R points. That determined the R = value of the dual fiuxoid in eqn.(p!9|), 
giving a prediction of e = 1.07(7), compared to the expected value 1/v^ = 1-03. All these 
features correspond directly to the classical solution of the lattice Higgs model reinterpreted 
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as a dual theory. 

Our goal here is to learn how to measure parameters of the dual superconducting medium. 
The Wilson loop projects out a normal region giving a normal-dual superconducting bound- 
ary and a chance to discover a London relation in the spatially transient region on the dual 
superconducting side. The U(l) model has a rather small correlation length for (3^1 and 
it is difficult to probe large distances. One small improvement is to replace the Wilson 
loop projector by a plaquette and study the tails of the correlators, Eqn(pO|). We present 
some preliminary observations but this small scale simulation is inadequate to tie down a 
consistent determination of A^. 

Fig.0 shows such a plot. This represents 2000 measurements for f3 = 1.0 and for a 
plaquette projector. An approximately exponential fall off is expected. The curve is only to 
guide the eye, it is an eyeball fit to the tail of the electric field: 

Q-R/Ro 

EiR) = Eo^^ (21) 

where Eo = 0.2 and i?o = 0.7. We use the 4-dimensional behavior for the tail for these 
'point like' operators. 

One can read off the point by point determination of the dual London penetration depth 
from Fig.0 



" V 27r[cMr/J-] ^^^^ 



On axis, at R = we have instead 



^ \l 2n[curlJ--E/{E-^^)] ^^^^ 

and the quantity in square brackets is plotted instead of [curlJ"^{0)] in order that the R = 
determination of can be read from the graph the same as the R cases. 

If we look along a lattice axis transverse to the projector plaquette , i.e. integer R, then 
E/curlJ"^ ^ 1.0. giving A^ ~ 0.4. compared to the eyeball fit of 0.7. Off the lattice axis, R 
noninteger, E/curlJ'^ ^ 2.0 giving A^ ~ 0.56. We expect the point by point determinations 
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of Xd should be consistent with each other and with the exponential tail determination only 
at large distances. This small scale simulation is inadequate to make a thorough check of 
this. 

A run of 2420 sweeps at /? = 0.95 gave similar results. An eyeball fit to the E data gave 
Eq = 0.4 and Rq = 0.5 which is consistent with the Wilson loop determination above. This 
agrees with the expected result that decreases with decreasing (3. 



IV. DUALITY TRANSFORMATIONS IN U(l) 

We arrived at the London relation in the U(l) theory purely numerically in Section 
|T|. This begs the question of whether the dual formulation of U(l) lattice gauge theory is 
equivalent to a lattice Higgs theory with a spontaneously broken gauge symmetry. Frolich 
and Marchetti have found this to be the case which we review briefly in this section. 



The duality transformation has not been performed for the Wilson form of the action 
which we are using. But this transformation has been implemented in closed form in the 



closely related Villain [20| form of the action by Frolich and Marchetti |23] and developed 



further by PoUey and Wiese, and Polikarpov and others i2^ , |25|j26| , |27|| . We draw on these 
later papers to gain more insight into the dual London relations. 

Recent literature on this subject have found it very convenient to use the differential 



form notation adapted to the lattice. PoUey and Wiese P4| , Sec. 2.1] have given a brief 
introduction to the subject. The advantage of this formulation is that it forces one to focus 
on the k-volume cell on the lattice for objects with k indices. Further it dispenses with the 
indices and reduces the typical manipulations to some very simple general operator relations. 
As a supplement to PoUey and Wiese's introduction, we give the reader a concrete realization 
of the differential form algebra in the appendix |A 1| in order to make the transition to this 
notation a little easier. 

Consider the Wilson action 
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Z = 11 r dip exp{-/3^(l - cosdifi)}, (24) 
■^-■^ dip 

where the hnk angle is a 1-form and the plaquette dip is a 2-form constructed with the 
exterior differential, d. We can expand each plaquette in a fourier series 

exp{-/5(l - cosz)} = J2e'''e~''h{P), (25) 

k 

^Y^^kz^-k-im^ (26) 

where we have used the asymptotic form of the Bessel function /fc(/?) for large /3 and further 



for large k. This gives the Villain |^ form of the action 

Z = n /_/^ HE expO(fc,rf^) - ^ P lin, (27) 
where the inner product, ($, \E'), is defined in the appendix |A 1| and || $ |p= ($, $). We now 



use the property of the inner product, {k,dip) = {6k, ip). We can now integrate over the 
links (fi and we obtain 

^ = nE exp{-i- II A: f } . (28) 

^ 2 (J 6k=0 

The variable k is an integer valued 2-form analogous to the field tensor F^^^, and the constraint 
6k = is analogous to d^F^^ = = 0. Next, using the general operator relation 5^ = we 
solve the constraint by introducing 

k = 6p, (29) 

^ = HE exp{-^ II 6p f } = nE exp{-^ II d*p f }. (30) 

The variable p is a 3-form analogous to a dual potential written as an antisymmetric three 
index object, Ba/s-y and Eqn. (|29|) is analogous to F^^, = daB^^y. If we switch to the dual 
objects this last result looks much more familiar: *p is a 1-form, *k = d*p *F^y = 

d;B, - diB,. 

In order to verify this interpretation we should introduce external sources and observe 
how they couple to the fields. Equation (^) naively looks like a free field theory, except 
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that the variables are integer valued and this indeed makes it an interacting theory. One 
can trace the integer value to the fact that the link angles are compact and hence their 
Fourier transform is a discrete variable. But this compactness is also responsible for the 
existence of monopoles. 

Our immediate concern in this paper is to show that this interacting theory given by 
Eqn. (|30|) is in fact a dual lattice Higgs model ^Sf. It is interesting how this comes about 



and gives us some insight into the significance of an integer valued dual gauge field. 

Following references pSj^ we introduce a non-compact 1-form gauge field *A on the 
dual links. We choose the Higgs field to be compact 0-form *\E' = exp{i*x} on the dual sites 
constrained on the unit circle. Finally we introduce an integer valued 0-form *i on the dual 
sites 

Z = l[r d*A UT d*x l[Y.^M-7^\\d*Ar-^\\d*x + 27r*i-*Ar}. (31) 

This differs from the lattice Higgs action, Eqn.(^^, in two respects: (i) The gauge field is 
non-compact, and (ii) the Villain PD| form is employed. The variables *i are introduced in 
order to make the action periodic in *x- 

To connect this theory to Eqn.(pil|) we first choose a gauge for which *x = 0. Then note 
that for large the integral over *A is peaked at the values *A = 2TT*i. In the limit k^i oo 
this becomes 

^ = nE exp{-|^ f }. (32) 

Identifying *i with *p this is the same as Eqn.(0) with 1//5 = = (27r)^/e^. 

Following Sec.|ITT| the dual London penetration depth in this theory is given by A^/ = 
l/(m^a)^ = l/(e^fi;rf). Therefore in the weak magnetic coupling limit, 0, the U(l) 

gauge theory with the Villain form of the action becomes a dual Higgs theory with an 
infinite photon mass, i.e. zero London penetration depth, A^. We find in our simulations 
that Xd decreases for decreasing (3. But we have only measured it for a very limited range 
of 0.9 <p < 1.0. 
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The dual Higgs field, *\& has been connected to the monopole field operator p3| , p4 . 



Polikarpov, Policy and Wiese have calculated the constraint effective potential for this 
operator and showed that there is a global symmetry breaking in the dual superconducting 
phase. 

V. GENERALIZATION TO PURE SU(2)AND SU(3) GAUGE THEORIES 



In ref. ||ri| we applied these same techniques and showed that dual Abrikosov vortices 
also occur in SU(2) pure gauge theory in the maximal Abelian gauge. Matsubara et.al. fl^ 
confirmed these results with better statistics and generalized them to SU(3). The SU(2) 
link matrices are U^{x) and the action is 

^ = /3 E {l-\TrU,,{x)). (33) 
The maximal Abelian gauge is defined by maximizing the quantity 

R = Y.Tr[a;U^{x)a^Ul{x)]- (34) 

The Abelian link angle is then taken as the phase of [f/^(x)]ii and the calculation can proceed 
with little change [^,0 . 



Recently we have done the calculation at finite temperature in order to check this picture 
on each side of the deconfining phase transition [l^ which we report here. The results are 
shown in the confining phase, fig.|^, (3 = 2.28 and the deconfining phase, fig.p, (3 = 2.40 on a 
lattice 4 x 17^ x 19, with 800 measurements for each case. Gauge fixing required about 600 
sweeps for each configuration. 

Fig.| shows an important difference from the U(l) case. Whereas in the U(l) case a 
linear combination of E and curlJ"^ for r 7^ can vanish giving the London relation, it is 
clearly not possible here. The behavior of —curlJ"^ does not match that of E. We interpret 
this discrepancy as a signal of a non-zero Ginzburg-Landau coherence length, ^d- Unlike the 
U(l) case, the extreme type II limit, in which the superconducting order parameter turns 
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on at the surface, the evidence here is that the order parameter turns on over a distance ^d- 
The value of the coherence length is approximately the radius where the London relation is 
restored. 

The analysis in this more general situation precludes a point by point comparison of the 
E and curlJ"^ data. We adapt the analysis from Tinkham [Q]. 

In terms of a Higgs model, the Higgs field is not constrained to a certain vacuum expec- 
tation value but rather is subjected to a SSB potential Eqn. (|T7|) . The equivalent GL theory 
provides a way to calculate the deformation of the Higgs field between the normal 
core of a vortex and the asymptotic value at large distance from the vortex, 0oo- We write 
= (f)oof{r)expia{x), where f{R) rises from — 1 over a distance ^d- Introducing a 
dual vector potential. Ad, the dual London equation is replaced by the dual GL equations 
for the azimuthal component of the monopole current. 

E = Vx Ad, 

J- = ^iAd--^Va), (35) 

and / satisfies 

-dVV + ei^a - fldYf -f + f = 0. (36) 
The generalized fluxoid relation becomes 

E{r) - A^Vx (:^) = N^J{x^). (37) 
An approximate solution of Eqn. (^5]) is 

f{r)=tanh{ur/^d), (38) 

with z/ ^ 1. 

An interpolation of the data is required to fit Eqn.(|37D. An example of the fit for j3 = 2.4 
in Ref. ^ is \d/a = 1.05(12) and ^d = 1.35(11). 

The behavior of all the SU(2) and SU(3) examples in the confined phase is similar to 
that shown in fig.^]. The interesting conclusion is that 
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t^d = \d/id ~ 1. (39) 

(no relation to Higgs k) For a type II dual superconductor Kd> 1/ \/2 and a type I otherwise. 
These simulations indicate that the non-Abelian dual superconductors lie at the borderline 
between type I and type II. 

Matsubara, Ejiri and Suzuki |T^ have measured for SU(2) for 3x3 and 5x5 loops in 
the range 2.4 < [3 < 2.6; and for SU(3) for 3 x 3 loop in the range 5.6 < (3 < 5.9. Although 
they are not asymptotic, they find clear signals in these data for both larger and smaller 
than 1/^2. 

The simulation [l^ shown in fig.^ was identical to fig.| except that (3 was increased to 



2.40 on a lattice 17^ x 19 x 4 putting us in the deconfining phase. The dramatic decrease of 
curlJ"^ results in the failure of the dual superconductor interpretation as expected for the 
deconfining phase. Also E falls more slowly with radius. 

A further observation is that although this dual superconductivity signal drops rapidly 
above the deconfining temperature, the monopole density does not change dramatically. The 
monopole density p is defined 

where \JTiat\ is the integer valued lattice quantities and Nsu^s is the number of lattice sites. 



In Ref. 1 13 1 we find 



confining : /? = 2.28 p = 0.0567(1) 

deconfining : /3 = 2.40 p = 0.0213(1) 



VI. EXTERNAL ELECTRIC FIELD 

Many interesting superconducting properties can be elucidated by studying the proper- 
ties of the material in the presence of a magnetic field. We report a preliminary look at 
the corresponding problem of dual superconductivity in the presence of an external electric 
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field. The Wilson loop provides such a field by projecting a qq out of vacuum configurations. 
But it may be interesting to see the spontaneous breaking of translation invariance as the 
electric flux forms a vortex rather than impose the vortex position at a given location and 
with quantized flux. Type 1 and type 11 dual superconductors respond very differently to 
background uniform external flelds. 

For a periodic U(l) lattice the sum of the plaquette angles over any plane is identically 
zero. Referring to Fig.[l^ a classical uniform electric fleld for U(l) on a particular time slice 
can still be obtained by constraining one plaquette angle in each z, t plane to the value 
{l/{NzNt) — l)Oc- Then the remaining {NzNt — 1) plaquette angles will take the common 
value 9c/ (N^Nt) = ea^F^^ giving a uniform electric fleld on all but one time slice. (N labels 
the lattice dimensions.) Choosing 6c = 7^ gives the largest fleld. For an 8^ lattice and 
/3 = 1 for example, there would be enough electric flux to form about 3 vortices in the 
X, y plane in the dual superconducting phase. This fleld conflguration can also be obtained 
by multiplying the plaquette that was singled out in the action by a minus sign and such 
conflgurations have been studied in non-Abelian theories |^^. We avoided an alternative 



method of imposing an external fleld by introducing a non-zero equilibrium value locally for 
each plaquette angle since we eventually want to see the fleld break translation invariance. 
Turning on interactions for (3 < 1 brings up other interesting features . 



Our goal here is to try to see a signal showing that curlJ'"^ responds to the external fleld. 
The immediate problem is that the sum of {curlJ^)xy over the any x,y plane is identically 
zero, and the sum of E^, is not. A sample conflguration is given in Fig.|lT[ Yet we expect the 
London relation to be satisfled. The only possibility is that translation invariance is broken 
which is expected since vortices segregate the superconducting and normal phases. 

We make the following rough ansatz that the local London relation is due to a the 
alignment of local current loops in the external fleld. This suggests that we truncate curlJ'^ 
to include only values ±2, ±3, ±4 representing the winding of the current around the dual 
plaquette. We denote this J™(2 + 3 + 4). In this discussion we assume that the current 
takes only the values —1,0, 1. The fractional population with values outside this range is 
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small enough to be neglected. 

Fig.^ is a typical configuration for the confined phase. A value ±4 indicates that the 
current makes a complete loop, the value ±3 a partial loop, etc. The value ±1 corresponds 
to a single isolated straight current segment in the plane. One expects the local behavior 
of ±4, ±3, ±2 contributions to align themselves in the local external field. But there is no 
comparable local behavior for the ±1 contributions. Globally of course they must all sum 
to zero. 

On a large lattice where size effects are irrelevant, one can still expect curlJm to average 



to zero. This is precisely the behavior in an isolated vortex as discussed in Sec. II B 



Now we get a large signal for curlJJ^{2 + 3 + 4) as shown in fig.|T2| with a sign that agrees 
with the sign in fig.§ for the superconducting region, r ^ 0. Further if we recalculate fig.|^ 
with the truncated current, we find a large suppression at r = and a moderate enhancement 



at the other points as shown in Fig.|T3|. In other words this choice biases in favor of the dual 
superconducting phase. 
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APPENDIX A: NOTATION CONVENTIONS 
1. Lattice differential forms 



The language of differential forms on the lattice is particularly useful for the Villain |2^ 



form of the action. Frolich and Marchetti p3i use this language in the mathematical physics 



literature. PoUey and Wiese Sec. 2.1] have given a nice introduction to the subject. This 
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appendix is meant only as a supplement to that introduction by relating the differential form 
operators to ordinary lattice finite difference operators. The differential form operators have 



coordinate independent definitions . The correspondences in this appendix allows one to 
check the algebra of differential forms on a specific coordinate system naturally defined by 
the hypercubic lattice. 

The lattice forward and backward difference operators are defined 

A+^{x) = ^{x + e^''^)-^x), (Al) 
A-^x) = ^x) -^x -e^^^). (A2) 

We offer the following examples of the use of the exterior differential d, and the codifferential 
6 acting on objects defined on the original (rather than the dual) lattice. 

0- form defined on points $(co) <^==^ <l>(x): 

^ A^^x), (A3) 
5^ 0. (A4) 

1- form defined on links $(ci) <^==^ ^f^{x): 

^ A+$,(x)-A+<l'^(x), (A5) 
5$ ^ A^$^(x). (A6) 

2- form defined on plaquettes $(02) <^=^ ^fiu{x) (antisymmetric): 

d^ ^ A+<|.h(x)= (A7) 

A+$^,(x) - A+<l>,,(x) - A+<l>^„(x), 
5$ ^ A^$^.(x). (A8) 

3- form defined on 3- volume $(03) <^=^ ^fii/\{x) (completely antisymmetric): 

d^ ^ A+$^,A](x)= (A9) 
A+<l>^^A(a;) - A+$«^A(a;) - A+<l>^«A(a;) - A|<l>^^„(x), 
^ A;^^,x{x). (AlO) 
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4-form defined on 4- volume $(04) <^=^ ^^u\t{x) (completely antisymmetric): 

d(!> 0, (All) 

5$ ^ A^<l>^,Ar(x). (A12) 

Miscellaneous relations on k-forms: 

52 = 0, (A13) 

= 0, (A14) 

Sd + d6 = A ^ (^15) 

<^ = {6A-^d + dA-^6)<^, (A16) 

(<l>,xl>) ^ 5: $,,,...(x)vl/^,...(x), (A17) 

X,fJ,>U-- 

($,c/^) = (5$,^). (A18) 

Note that the two arguments of the inner product must of course live on the same k-form. 

There is a completely equivalent description of these k-forms, $(cfc), in terms of (4-k)- 
forms on the dual lattice *$(*Cfc). 

*$(*c,) = <l>{ck). (A19) 

This is realized in conventional notation by contracting e^i^xr into the object defined on the 
original lattice to obtain the object on the dual lattice, properly normalized so that acting 
twice is the identity operation. 

Finally we wish to note a correspondence involving d and 6. Given the equivalence 
Eqn. (|A19|) then there is a correspondence between the {k ± l)-form and dual (4 — (A; ± 1))- 
form 

d^ = 5*$, 

5$ = (A20) 

On the original lattice the position x refers to the 'lower left etc' corner of a k-form. 
(Links point out of their coordinate address.) This dictates the choice of A"*" and A~ in 
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these examples. Since the dual lattice is displaced a half spacing in the positive direction, 
the 'upper right etc' corner is used to label the position x of the dual k- forms. (Dual links 
point in to their coordinate address.) As a consequence, when the above rules are applied 
to k-forms on the dual lattice, A+ and A~ should be interchanged. 



2. Euclidean Maxwell equations 

We have taken the following sign conventions for the Euclidean field tensors with electric 
and magnetic sources. 



d F — 




* Z? 
^ap — 


— -^^aPixvFn^, 61234 — 1, d*^FiJ.v — 






V X B- 


- diE = fi, 




jm 
^4 , 


V X E 


- diB = J^. 



For — and with the standard definition F^^ — d^A^ — d^A/j,, we get the unconventional 

— * — * — * 

sign B — —V X A. 
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FIGURES 



FIG. 1. DeGrand-Toussaint construction to identify a monopole in a spacial 3-voIume. 

FIG. 2. (a) Six space-time plaquettes (time not shown) contiguous to a time-like link, (b) Six 
space-space plaquettes dual to those in (a) forming a cube. 

FIG. 3. Profile of ea'^Fuix), Eqn.(|l|), and -(aVeK)(A+ J|(x) - A+Jf(x)), Eqn.(||), as a 
function of distance perpendicular to the plane of the monopole current loop in lattice units. 



FIG. 4. Operators for (a) ea^Fsiix), and (b) (aVem)(A^ J2™(x) - A^Jf (x)). 

FIG. 5. Profile of ea^F34(x) and — (a^/em)(A|f J2"(x) - A2"J{"(2;)) as a function of distance 
perpendicular to the plane of the Wilson loop in lattice units for U(l). 

FIG. 6. Profile of ea^F34(x) — A^(a^/em)(Aj~ J2"(x) — A2"J{"(x)) as a function of distance 
perpendicular to the plane of the Wilson loop projector in lattice units for U(l). 

FIG. 7. Profile of ea^-F34(x) and -|-(a^/em)(Aj~ J2"(x) — Ag J{"(x)) as a function of distance per- 
pendicular to the plane of the plaquette projector in lattice units. The on-axis point for cur/ J™ < 
is 'corrected', see text. 

FIG. 8. Profile of ea^F34(x) and — (a^/em)(A||" J2"(x) — A^J{"(x)) as a function of distance 
perpendicular to the plane of the Polyakov lines separated by i? = 3 in lattice units for SU(2) in 
the confined phase. 

FIG. 9. Profile of ea^F34(x) and — (a^/em,)(A^ J2"(x) — A^J{"(x)) as a function of distance 
perpendicular to the plane of the Polyakov lines separated by i? = 3 in lattice units for SU(2) in 
the unconfined phase. 



FIG. 10. Plaquette angles for a classical external field configuration. 
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FIG. 11. Sample values of the z component of curUJ^^ on an x,y plane for fixed z and t on a 
lattice, /? = 0.99. 

FIG. 12. Average of 60^^34 and -(aVem)(Ai J^(2 + 3 + 4) - A2 Jr(2 + 3 + 4)) on each time 
slice for an external field = 7r/64 (the horizontal line). The constrained plaquette is at t = 1, the 
field is classical at t = 2, 8 and /3 = .99 for t = 3 — 7. 

FIG. 13. Profile of -{a^ /em){^i J^ix) - Jr(2;))(2 + 3 + 4) as a function of distance per- 
pendicular to the plane of the Wilson loop in lattice units for U(l) superimposed on Fig.^. 
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